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SUMMARY 

A new explicit method for solving the transonic sma*!- 
disturbance potential equation is presented. This algorithm, 
which is suitable for the new vector-processor computers such as 
the CDC STAR-100, is compared to successive line over-relaxation 
(SLOE) on a simple test problem. The convergence rate of the 
explicit scheme is slower than that of SLOR. However, the 
efficiency of the explicit scheme on the STAR-100 computer is 
sufficient to overcome the slower convergence rate and allow an 
overali speedup compared to SLOR on the CYBER 175 computer. 

INTRODUCTION 

The State-of-the-art of transonic flow calculations has 
advanced to the point where two-dimensional flows, including the 
effects of viscosity, can be computed in a relatively short time 
on modern serial-type computers. For example*, many people are 
using a program developed at the Couraiit Institute of New York 


if. 

Courant Institute of Mathematical Sciences, New York University 



University for the analysis of transonic flow past airfoils 
(ref. 1). This program gives accurate solutions to the full- 
potential equation, including the effects of boundary-layer 
displacement, in about 2 or 3 minutes on a CDC CYBER 175 computer. 
Three-dimensional, transonic, finite-difference calculations, 
however, are expensive on this type of computer. The three- 
dimensional program described in reference 1 takes about half an 
hour for an inviscld calculation on a fairly crude grid. 

It is hoped that the use of the STAR-100 computer will allf)w 
accurate, three-dimensional, transonic flow calculations to be 
done economically. One way to achieve this goal is through the 
development of algorithms which can make full use of the unique 
architecture of the STAR-100. The STAR computer has a "pipeline" 
type of processor which is very efficient in doing arithmetic 
operations on long vectors (ref. 2). Unfortunately, the best 
available method for solving the transonic potential equation is 
successive line over-relaxation (SLOR), which is not amenable to 
vector arithmetic. The reason for this is the semi-implicit 
nature of the iterative method; that is, the calculations at a 
particular grid point require results from the current iteration 
at neighboring grid points and thus cannot be done in long vector 
operations. 

This paper describ'es a new explicit algorithm which can be 
vectorized for use on the STAR-100. The new algorithm is applied 
to a simple test case and compared to SLOR on the CYBER 175 


computer. A simple STAR program, including vector instructions, 
is given in the appendix. 


TEST PROBLEM 

The test problem chosen for this preliminary study is to 
solve the transonic, nonlinear, small-disturbance potential equa- 
tion for a nonlilting parabolic-arc airfoil in a finite box with 
uniform grid, as shown in the sketch below. 
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Although this is a very simple physical situation, it still has 
some of the most difficult features of transonic flow fields as far 
as programing for the STAR is concerned. The governing partial 
differential equation is 

1 - - (y+ 1) d) Id) + d) =0 

L 00 00 't'xJ 'xx ^yy 

The boundary conditions to be applied are 
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{(> = 0 on the outer boundary 


and 

(J) = +4tx on y = ±0, -.5 ^ x ^ .5 

i/ 

where t is the thickness of the parabolic-arc airfcil (which has 
a unit chord). For t = 0.1 and = 0.9 the flow is super- 
critical. In regions where the coefficient of (j> is positive, 
the flow is subsonic and the equation is elliptic type. In regions 
where the coefficient of (j> is negative, the flow is supersonic 
and the equation is hyperbolic type. The general procedure for 
solving this equation is to replace the partial differential equa- 
tion with a finite difference equation at each grid point. These 
finite difference equations are then solved iteratively. 

This test problem represents a simple physical situation which 
is of little practical interest. A more useful program should 
allow for lifting flows and should extend the outer boundary 
farther away from the airfoil. This could be done either by using 
a stretched grid or by some type of grid nesting using additional 
coarse grids around the small region considered here. The test 
case does, however, include the major difficulties to be overcome 
in using the STAR-100 computer for transonic flows. For example, 
it has supercritical flow which requires a change from one type of 
difference equation at subsonic (elliptic) points to another type 
of difference equation at supersonic (hyperbolic) points. It also 
requires the use of an explicit Iterative scheme to solve the 



difference equations if the arithmetic operations are to be done 
using long vector Instructions. 

SEMI-IMPLICIT SOLUTION METHOD 

The most common method used to solve the finite difference 
equations is successive line over-relaxation (SLOP). This itera- 
tive scheme is implemented as follows: 

Compute U 5 1 - M^ - (y+1) (|>^ 

ra uQ ^ 

using 


n 




X 


= ^‘i-H.-l “ ^i-1.,1 
2Ax 


If U > 0 (subsonic points), central differences are used 
together with over-relaxation (w > 1) to give: 


U 


.n 2 .n+1 

'l’i+l,d " 0) '^i,J 



1 

(I) 





n+1 

i“l,j 


,n+l 

i..1+l 


- 24.V 


n+1 


n+1 




lij-i 


Ay' 


0 


where the superscript indicates the iteration number. 

If U < 0 (supersonic points), an upwind difference is used 


for 4 to give: 






At the airfoil boundary tho 4*^ is replaced by 

^ J 



which takes into account the fact that tho nonlifting flow is 
symmetric (i.o., at y = 0). Each iteration is 

generated one column at a time going from left to right, by 
solving a tridiagonal system of equations at each column. This 
nonconservative scheme is similar to that originally proposed by 
Murman and Cole (ref. 3). See references 4 and 5 for a discussion 
of related conservative schemes. 

EXPLICIT SOLUTION METHOD 

The explicit solution method uses values of the potential 
function from the two previous iterations in ordei’ to update the 
potential function for the present iteration (thus termed a three- 
level scheme). The method is implemented as follows: 

Compute U = 1 - - (y+l) 0 

CO CJ 

using 
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If U > 0 (subsonic points), central differences are used at 


iteration, or "time," level n: 




n+1 


- ‘’■♦i.j *5:5 ^ (2 - P] 



- <!> 


i,J/ 




where 


j - 


j 


i-l, j 


M! 

Ay" 




n 


iJ+1 


“ 2i\> 


n 


i) j 




and 


D 


1 



If tf < 0 (supersonic points), the upwind <j) difference 
is formed at "time" level n - 1, while (p^ is centrally~ 

•7 «/ 

differenced and averaged in time level n, to give*. 


.n+1 ,j.n . .n-1 . 

"^i.J “ ^^i,J ^i,J ^°\'^i,d " *^i,j ^ ^i“l,j ^i-l.j 


lU-1 


n~l 


, n-1 


“ ^ 2^2 


( 1 ) 


Where 
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^2 = u j) 

* “(♦i-l.j+l ■ ^'•'i-lj * 'I’l-l.j-l)] 


U “ (1-a) U + 0U 


U - 1 - - (Y+1) 


,n-l ,n-l \ 
^1..1 - 1’l-2..i ) 
2Ax / 


0 = Pg min 



and 



A von Neumann stability analysis of this last scheme with 
the U replaced by U shows why the tj) derivative is evaluated 

aA 

' ' w 
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at the n-1 iteration and why the i|)yy derivative is a weighted 
average betv/een (|)yy at the i column and <})yy at the i~l 
column. To do the von Neumann an;V*.ysis, let i|) at the kth 
iteration be 

Also let mAx = C nnd nAy = n. Putting these definitions into 
equation (1) gives 

g^ - 2g + 1 + 2o (g - 1 - ge”^^ + e"^^) 


= ^2\ 

U (1 - 2e-« + 


t e 

(Hf [a-o)(e^'’ - 2 + e'^") 


+ a 


(2) 


Note that 


DgU = - 


Also define 


l-OT (#) 5 
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so equation (2) becomes 


- 2pg [x - a (1 ~ e"^^)] + [l - a (1 - ■= 0 

or 

g “ (p i i six - P^)[l - a (1 - e"^^o| 

Thus, in order to have |gj - 1. it is necessary and sufficient tn 
have both |p| ^ 1 and a £ X. The inequaXity |p| ^ X implies 
that a sVTuf 

This convenient factoring of the expression for the amplifica- 
tion factor was made possible by choosing this particular repre- 
sentation for the and <Jiyy derivatives. In practice the 

coefficient of (fi is tf rather than U; this averaging makes 
the scheme approach second-order accuracy as a approaches one. 

It should be noted that this explicit scheme is not only a 
different iterative algorithm from the semi-implicit scheme, but 
it also has a different steady-state solution. This difference 
occurs because of the weighted averaging done on both the <j> 
term and on the coefficient of the <j)„„ term for supersonic 
points. 

CONVERGENCE AND TIMING COMPARISONS 

Short computer programs have been written to solve the sample 
problem using the two methods described. The semi-implicit SLOR 
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method was coded in standard FORTRAN IV and run on a GDC CYBER 1,75 

1 

computer (which is about 2 ^ times as fast as a GDC 6600 for 

this type of problem). The explicit method was coded in STAR 

FORTRAN (ref. 6) and run on a GDC STAR-100 computer. The STAR 

code Includes the use of vector instructions in the iteration 

loop. It also includes the use of bit control vectors to distin- 

(juish between subsonic and supersonic points. The bit control 

vectors provide the capability of performing the complicated 

supersonic calculations only at supersonic points, which are 

collected into a vector through the "compress'* and "expand" type 

of instructions available on the STAR. The STAR code for this 

problem is listed in the appendix. 

The computations wore done on three different grids. Each 

calculation was terminated when the value of the lai’gest residual 

12 2 

in the flow field was less than g + Ay ) , Each calculation 

was run with experimentally-determined optimum values of the 
parameters for that algoritlim so that convergence was attained in 
a minimum number of cycles. The results in the table below show 
that the new three-level explicit scheme has a slower convergence 


Grid 

size 

SLOR 

ON CYBER 175 


EXPLICIT METHOD ON STAR-100 

Cycles to 
converge 

Time to 
converge 

Average 

sec/cy. 


Cycles Uo 
converge 

Time to 
converge 

Average 

sec/cy. 

40X40 

42 

1.242 

.0296 

n 

131 

.879 

.0067 

80X80 

98 

12.118 

.1237 

1 

300 

6.385 

.0213 

160X160 

244 

116.095 

1 

.4758 j 

1 

u 

1 655 

i 

59.834 

,0913 
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rate than SLOU. IIowe»vor, the offlcioncy of this now schomo on tho 
STAR computer is enough to make up for the slower oonvergenco rate 
and still allow an overall reduction in computing time. Not 
surprisingly, the speedup is greater for tho cases with finer 
grids than for the 40X40 case, 

CONCLUDING REMARKS 

This preliminary study has shown that a new explicit method 
for solving the transonic small-disturbance potential equation on 
the STAR-100 computer can almost halve the computer time required 
for this type of computation when compared to successive line over- 
rolaxation or- the CDC CYBER 175 computer. These results are 
limited tc u. relatively simple problem with a uniform cartesian 
arid. Although the speedup is not as great as desired, it is 
enough to justify further study of this method. The effects of 
lift and of grid stretching on the convergence rates of the 
schemes should be investigated. Also, the new explicit scheme 
(whould be applied to the full potential equation. 

There are several possibilities for obtaining further reduc- 
tions in computer time. One is through the development of more 
efficient algorithms. Improvements might be made in the conver- 
gence rate of explicit algorithms, or other vectorlzable algo- 
rithms might be developed. Another possibility is through 
programing techniques to get successive line over-relaxation to 
run as efficiently as possible on the STAR-100 computer. Although 
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the method is semi- implicit , there are portions of it which can 
bo written in vector instructions of short length. 

The program listing in the appendix should serve as an 
introduction to some of the programing techniques available on 
the STAR-100 which are useful for transonic flow calculations. 


APPENDIX 


PROGRAM LISTING 


PROGRAM UT ( I f«PUT» OUTPUT /TAPES* I MPUT/TAPE6aOltTPUT) 

COMMOK* P(/il ,PP{/a ),IMAX,JMAXf J»«/CnVRrMiT,lMAXMl,Jf'AXMl 

* ,nX,OY,OX?R;Pl ;P2,X*^lNF,RSH,DYaR 
CALL R3CLOr.KS(CPil/^«ALL) 

KKnE(6,P9(>2)CPU,WALL 

XMjNFe,9 

Pl=t ,03 

PP.c.99 

PSHsl.a 

IHAXBill 

JMaXb/JI 

MITsl 000 

CFel. 

Pa=0,/3, 

PMsO,/3, 

IMAXMlsIMAX-1 

JMAXM1sJ.mAX«1 

DXsPK/IMAXHi 

OYsPH/JmAXmj 

JKtsjMAXMj/2 + 1 

0X2RKt,/DX*A2 

DYSRbI ,/DY**2 

C0VRcCF*,5*(0X**2+0Y**2) 

DO JO I«1/IHAX 
FP(I)aO, 

JO CONTINUE 
XLEs-,5 
XTEa,5 
EPSsl ,0E-06 
lHsINAXMJ/2-M 
ILEelFIXt (XLE“EPS)/DX)+IH 
nEelFIX((XTE+tPS)/DX5+lH 
DO 11 IkILE/ITE 
X=a"IH)*DX 
FPcnft<,o*x 
11 CONTINUE 

NTaIMAX*JHAX 
PU,l|NT) = 0, 

CALL PICTURE 
CALL PRESS 
STOP 

9902 FORMATdX, AlO, IX, AlO) 

END 


original page is 
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SUhRDUTINE PICTURE 

C THIS SUBROUTINE UOFS THE RELAXATION 

COMMON PC^Jt ,«1 )#FP(iin,lMAX,jMAX,Jw,COVR,MIT,lKAXHl ,JSAXhl 
* ,DX,0YfDX2R#Pl ,P2,XrlINF,R5H^pY2H 

C DIMENSION SOf'E ARRAYS FOR TEMPORARY STOWARE 
0 1 f'EMS 1 ON PP ( a 1 » /J 1 ) # pH ( « 1 , /J 1) 

DIMENSION TtCRl*^Jl)#T2(ai,«nrTj(/ll*«n 
BIT » IR(ai »NUD,Ef>D,IRO/BlD 

DESCRIPTOR T1 n, T?0,TJn/Pn,PPI)/PHD 
DESCRIPTOR Mun^EPD, IrtD^BlD 

DESCRIPTOR T lN'D,T2ND,T3Nn,T«ND^T5NP/T6Dn#T7ND#T6ND 
CALL 05CLOCKsS(CPU,i‘ ALL) 
write (6/RD0«)CPU,K all 

C ASSJRN THE OfSCRIPTORS TO SPECIFIC ARRAYS 
\>PeIHAT*JNAX-2*J"AX«2 
ASSIGN TlD^Tl (2,2>NP) 

ASSIGN T2D,T2(2,2M P) 

ASSIGN T3n,T3(2^?>f'P) 

ASSIGN Pn,PC2,2jNP) 

ASSIGN PPO,PPC?,?r'P) 

ASSIGN PND,Pm(S,?jLP) 

ASSIGN NUD,HU(2/2jNP) 

ASSIGN EHD,EIU2,2)NP) 

ASSIGN ibn,io(2,2ir.P) 

C SET UP A PIT ARRAY THAT HAS ONES AT INTERIOR POINTS 

C and ZEROS AT THE EDGE POUTS 
IRDnS'O' 

DO 1 Is2fU‘AXMl 
00 1 Jc2,JPAXMJ 
I IHCIrJ)=BM' 

EBDs.NOT^IPD 
XK=1 ,«XHINFA*2 

AC = (BSH + 1 , )*XMINF**2/(2,*DX) 

D1=0.5*PUP2*P2 , 

AR2SCOX/O.Y) A*2 
AR2RsJ ,/AR2 
02sP2*P2/AR^^ 

TW0DY52.ADY 

S=Pl-l, 

C INITIALIZE THE TEMPORARY ARRAYS 
TlDaO. 

T2DbO, 

T3Dsn, 

PPOsO, 

Pi'.DsO, 

V.RITE(6,9R0n) 

P0JNTSb(IM4X-2) A(J^‘AX«2) 

RNAXO=10, 

RAVGOslO, 

RASUB0=1 Oo 
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RASUPO»10, 

Np 1 

50 CONTINUE 

c compute U 

T10bXK-AC*(P(3#2jNP)-P(1,2»NP)) 

C SET THE Mun,J) BIT TO ONE AT SUPERSONIC POINTS 

HUDkT id.lt.o, 

C SET MU TO ZERO AT THE EDGE POINTS 
MUOsMUD.ANO, IRO 

C compute a CENTRAL PHi-XX AT IrJ 
PPOsPH ,2jNp)-PD-Pr+P(3#2|NP) 
c COMPUTE A CENTRAL PHI-YY AT IfJ 
TcDsP(g,l pNP)«Pn-P0+P(2i3>NP) 

C COMPUTE THE RESIDUAL 
PPDsT10APPD+AR2*T20 

C CHANGE the residual AT THE AIRFOIL SURFACE 

PP(2,JlMl^AX«2)BPP(2f J^.> I^AX-2 )+TVsDOY*FP(2j1MAX-2) 
T3P80. 

C PUT ZEROS IN THE RESIDUAL ARRAY AT THE EDGE POINTS 
PPDsR8VCTRL(T3D,EBD|PPD) 

C PUT ZEROS IN HE RFSIDUAL ARRAY AT SUPERSONIC POJNTS 
PPDsf;8VCTRL(T3D,MUDjPPD) 

T3DcVAHS(PPO|T3D) 

PSUPSUBsQ8SSUM(T3D)ADy2R 
RMAXSUB=nSSHAX (T3D3 *f>X2R 
C COMPUTE THE NEW PHI VALUES 

PPDs01/'(TlD + AR2)*PPD + PD-n ,"Pn*(PD-PMl>) 

C number of supersonic POINTS 
NS=Q8SCNT(MUD) 

RSUM.SUPsO, 

RMAXSUPaO. 

IF(NS,EO,0)GO TO 300 
C SUPERSONIC CALCULATIONS 

C ASSIGN temporary STORAGE FOR SUPERSONIC CALCULATIONS 
ASSIGN TlNDf.DYN.NS 
ASSIGN T2NT>» ,DYN,N8 
ASSIGN T3N0,,DYM,NS 
ASSIGN T«ND#,DYN,NS 
ASSIGN T5ND/,DYN,NS 
ASSIGN T6ND^ ,DYN,NS 
ASSIGN T7ND, .DYN.NS 
ASSIGN TSNDi.OYN.NS 
ASSIGN PID,,DYN,NS 

c make a vector out of the values of u at supersonic Points 

TlN0s06VCHPRS(TlD,Mun)TlNP) 

C COMPUTE SIGMA 

TSNDkVABSCT1mdjT 3ND)*AR2R 
TilNOcVSOkT CT3NDfTRND) 

BlDaT^JNn.GT, ! a 

TSND»1, 
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Ti!J»''()BP 2 *Q 8 VCTRL(T 5 f.*D,niD>TttND) 

C make a vector UUT np THE VALUES OF THE CtwTRAL 
C PhI"YY AY SUPERSONIC POINTS 

T2NDsOflvcMPRS(T2n,Mun»TaNO) 
c make a VECTOR OUT OF THE VALUES OF THE CENTRAL 
C Phl-YY ONE POINT IIPSTRE-M OF EACH SUPERSONIC POINT 
T5Nnsf30VCMRRSnP( J ,2iNP),MUDjTSN05 
C COMPUTE PHI-YY TO USE AT SUPtPSONTC POINTS 
T 5 NDsT 2^0 + T‘!^^■^*(T 5 ^ D-T 2 ND) 

C make vectors out of OLO values OF PMI AT SUPERSONIC POINTS 
c AND OWE AM) TWO POINTS UPSTREAM OF EACH SUPERSONIC POINT 
TSNDBOSVCMPkS (PMO,>‘U0f T6MD) 

T7(‘ObOPVCHPRS(PM{ J ,2f NP) ,MUI)>T 7Nf)I 
TKMDs()flVCMpRS(PM(0,2jNP) ,HUDiT8ND) 

C COMPUTE U«TILOE 

TgNpsXK-AC* (TSND-TflND) 

C COMPUTE U«bAR 

T2NDsTlNO-T-!JNO*(TlNn»T2Nn) 

C COMPUTE THE RESIDUAL 

T5NDsT5N0 + T2N0*n6ND“T7ND«*T7N0 + TSN|)) 

C EXPAND THE HESIOUAl.S AT SUPERSONIC POINTS lUCK 
C TO A TEMPORARY FULL ARRAY 

T30«0eVXPNnCT5NO,MuniT30) 

C change the RESIDUALS AT THE AIRFOIL SURFACE 

T3(2/ Jwf IMAX-2)=T3C2, JWj lMAX-2)+TWOPY*FP(2> IMAX-2) 

C REC0HPRES5 the. RESIDUALS AT SUPERSONIC POINTS 
T5MDn0aVCMPRS(T3P/M(jD|T5NO) 

C COMPUTE D2 

TlNO = VABS(T5WOrTU'D) 

PSUMSUPsOaSSUMCTl nd)*DX 2R 
RMAXSUP=NaSMAX(TlNO)*DX2R 
TlA'Osl , 

TlND«D2/UBVCTRLn3MD,BlD|TiND) 

C COMPUTE the new VALUES OF PHI AT SUPERSONIC POINTS 
TBND=0SVCMpR$(PDrMUD|T8ND) 

T2ND=T8M0-TbND 

T3NCs08VCHPRS(P(l ,2lNP) rMUD)T3ND) 

T3NDKT3ND-T7ND 

T5NDsTJNn*T5ND + T0NO + T2N'D«(T«ND + T<<ND)*(T2MD-T3NO3 
C EXPAND THE NEh VALUES OF PHI AT SUPERSONIC POINTS TO A FULL 
C ARRAY AND PUT THEM IN THE PHI-PLUS ARRAY 
T30 b08VXPND(T5ND,MUDjT3DI 
pposyevcTRLnsD^MUDrPPD) 

FREE 

300 continue 

C MOVE THE PHI ARRAYS TO N£W ITERATION LEVEL 
PMDsPO 
ppspPO 

PSUMsRSUM-SUB + RSUMSUP 
RHAXsAMAXI (RMAXSUB/RMAXSUP) 
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NTE»*PsNS 

IF (N5,ED,0)NTt»'P«l 
FNScFLOAT{K^TtHP) 

RASUPsRSUMSUP/FNS 

IF(NS.E«,0)WASHP3in, 

PASUBsCf{SIJM-.RSU^‘5UP)/(P0I^^lTS-FNSI 

RAVGsPSUf^/POlNTS 

SRt'sRk^AX/PMAXO 

SPAsRAVG/MAVUO 

SfiASimBPASUK/PASliRO 

SOASliPsRASUP/WASl'PO 

PMAXObRMAX 

P AVGOeRAVG 

RAStiHOBrnSUB 

RASUf OsRASUP 

i> PITE(F,9'?fla)N,RPAX,SP«f PAVG;5RA/RASl'n;SPASUR,PA5UP,SRASUP,NS 
lP(Ri^AX,LF»CDVR,nP,M,Ge,Mn )GO TO 500 
JF(^‘,F.Ra nPHAXIsK^'AX 
IF(^^tQ,nPAVGI=iPAV6 

Ni:W + i 

GO TO 50 
500 COnTIRIJE 

SR«a(R'^AX/K>''AXl)**n ./FLOAT CN)) 

Sk/ e (Pa VG/P AVGn**n. /FLOAT (N)) 

V.RITE(6,99('3)5PR/SRA 
CALL RSCLOCKSCCPI'^walL) 

VPlTE(6f9905)CPU,lvALL 

RETURN 

9901 FOR^'AKJHO N9XOHPMAX6X3HSRM5XOMRAVG6X3MSRA/JX5HRASUB5K 

* 6HSRA5UbaX5HR4SL'P5X5HSRSUP3X2HMS) 

9902 FORMATnX,I3,Ell,5,F7,9,E11.5#F7«0,El! ,5,F7.«rEll.5#F7,0, 10) 

9903 FORf-'AT('OOVEftALL SPECTRAL RADIUS I5'F7,0,' BASED ON RNAX AND' 

* F7.0#' BASED ON RAVG.O 

9900 F0RHAT(»0£NTERING PICTURE AF TER 'F 1 0 . 5 , ' SECONDS CPU TIME AND' 

* FIO.S,' SECONDS FLAPSED TIK£') 

9905 FOflHAT(' ITEPATinrts TDOK'FIO.5^' SECONDS CPU TIME AND'FlO.5, 

* ' SECONDS ELAPSED TIME') 

END 

SUBROUTINE PRESS 

common PC4WPl),FPC/in , IMAX, jHAX»Jt^;COVR,'^lT, IMAXMl r JHAXmI 

* #PX^DYfOX2R,Pi ,P2,XMINF,RSH,DY2R 
DIMENSION CP(ill) 

Js J W 

DO 1 Iog,IMAXMl 
CP(I)n(P(I«J , j)«P(1 + 1 ,J) )/0X 
I CONTINUE 

WR1TE(6,9902) (I ^CP(I), I=2,IMAXM1 ) 

RETURN 

9902 FORMAT('l I CP '/ ( I fl , F 1 0 . 6 ) ) 

END 
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